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I summarize recent developments in spectrum calculations using Kogut-Susskind quarks. Theoretical devel- 
opments include one-loop computations with improved actions. I present some recent simulation results, mostly 
from a MILC collaboration project using three flavors. Effects of dynamical quarks are clearly seen in the isovector 
0^"'' meson propagator and in the mass ratio "J" . 



The Kogut-Susskind formulation of lattice 
fermions, often called "staggered fermions", has 
long been a popular formulation for QCD sim- 
ulations including dynamical quarks. The basic 
reason for this is the single remaining exact chi- 
ral symmetry, which is sufficient to protect the 
quark mass from additive renormalization. From 
a practical viewpoint, an important consequence 
is a lower limit on the eigenvalues of M^M, which 
insures that the simulation will be well behaved 
and not encounter "exceptional configurations." 
Also, but less importantly, there are simply fewer 
fermionic degrees of freedom to handle. The com- 
bined effect of these advantages is that we can use 
our limited computer power to push to smaller 
dynamical quark masses, or larger physical sizes, 
or perhaps to larger statistics, by using Kogut- 
Susskind quarks. The price that we pay is that 
it can be a painful exercise to figure out the lat- 
tice implementation of your desired operator and, 
more importantly, flavor symmetry is broken by 
effects of order a^. 

Although they are of order a^, in practice the 
effects of flavor symmetry breaking are unpleas- 
antly large at accessible lattice spacings. The 
situation can be greatly improved by the use of 
an improved action which suppresses the coupling 
of the quarks to high momentum gluons, as well 
as fixing up the quarks' dispersion relation. See 
Ref. § for some relevant references. The most 
important point is that exchange of a gluon with 
momentum near n/a can scatter a low momen- 
tum quark from one corner of the Brillioun zone 
into another, resulting in a mixing of different fla- 
vors. "Fattening" the links, by averaging paths in 
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Figure 1. Mass of the quenched p in units of ri 
at TOtt'^i = 0.807. The squares are the improved 
Kogut-Susskind action[^ and octagons the con- 
ventional Kogut-Susskind action The plusses 
and fancy crosses are clover- Wilson action 0] and 
diamonds are domain wall quarks |8|. 



the parallel transport, effectively puts a form fac- 
tor into the quark-gluon vertex which suppresses 
such exchanges. 

At this conference we have seen the first re- 
sults of one-loop calculations done with improved 
Kogut-Susskind actions. In Ref. [|| the quark 
mass renormalization is calculated, and in Ref. |^] 
renormalization constants for several operators 
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Figure 2. The static quark potential with 
a ~ 0.13 fm. The octagons are the quenched 
potential and diamonds the three flavor potential 
at ms. The lines are fits to "Coulomb plus lin- 
ear plus constant" , and the rulers show the lattice 
spacings and units of 0.1 fm. 



are shown. The results are encouraging, with 
none of the coefficients coming out surprisingly 
large. In contrast, the conventional Kogut- 
Susskind action often has large one-loop correc- 
tions. This can again be understood in terms of a 
form factor suppressing coupling to high momen- 
tum gluons|^,^, so that a diagram like Fig. 2b of 
Ref. Q will not get large contributions when the 
quark is in another corner of the Brillioun zone. 
In other words, we suppress the unwanted dou- 
blers in the loop diagrams. 

Figure shows the result of one scaling test 
of the improved Kogut-Susskind action. In this 
figure we plot the p mass in units of the static 
quark potential as a function of lattice spac- 
ing. For the length scale we use ri, the distance 
where rfF{ri) = 1.0. All these points are from 
quenched simulations, and they are all interpo- 
lated to the quark mass where m^ri = 0.807, or 



Figures, rov^ as a function of quark mass. The 
quenched limit is at the right side. The octagons 
arc MILC results with three degenerate flavors at 
a « 0.13 fm, squares are two light quarks and one 
strange quark at a « 0.13 fm, and the diamond 
is a two flavor run at 0.13 fm|l^. The crosses are 
preliminary results at a « 0.09 fm. The bursts 
are CPPACS results with two flavors of Wilson 
quarks ||lj] and fancy diamonds are SESAM two- 
flavor Wilson results [p^. Some very new two- 
flavor Wilson results on smaller lattice (not plot- 
ted here) which show a much smaller effect can 
be found in Ref. ill. 



« 460 MeV. 
Next one is motivated to ask whether it is prac- 
tical to further improve the Kogut-Susskind ac- 
tion. At this conference DiPierro and Mackenzie 
reported on experiments with empirically tuning 
the coefficients of the paths in the action The 
gains were limited, which they interpret as evi- 
dence of the need for four quark operators. In 
Ref. Q a one-loop improved action was presented, 
and methods for handling these four quark oper- 
ators by introducing auxiliary bosonic fields were 
sketched. Taking a different approach, A. Hasen- 
fratz and Knechtli combined fattening of the links 
with a projection back onto unitary SU(3) matri- 
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Figure 4. Squared pion mass divided by quark 
mass. The octagons are quenched a « 0.13 fm 
results, diamonds two flavor, and squares three 
flavor. The bursts are preUminary quenched and 
three-flavor results at a w 0.09 fm. 
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Figure 5. Quenched Goldstone pion masses from 
Kim and Ohta0 and MILC||l|. The axes have 
been arbitrarily rescaled to match the values at 
the strange quark mass, which is approximately 
the low point on the curves. 



ces in the "HYP" , for "Hypercubic blocking" ac- 
tion[0. This action produces smaller mass split- 
tings among the pions than the tadpole action. 
These authors also describe an algorithm for dy- 
namical simulation of this action, which cannot 
be expressed in a simple way as a sum over paths. 

In another theoretical development presented 
at this conference, Levkova and Manke have 
worked out an action for unimproved dynami- 
cal Kogut-Susskind quarks on an anisotropic lat- 
tice While their primary motivation is high 
temperature QCD, this approach could be useful 
for spectroscopy of glueballs, hybrids and excited 
states. 

The MILC collaboration has used the tad- 
pole improved action for a set of hadron spec- 
trum calculations |lj]. These calculations used 
three flavors of dynamical quarks, as well as a 
quenched run and one two-flavor run, on lattices 
tuned to match the lattice spacing at about 0.13 
fermi. For masses larger than the strange quark 
mass nis three degenerate flavors were used, while 



for smaller masses one quark mass was held flxed 
at TTis. Quark masses down to mq = 0.14mc,, or 
mTj/nip — 0.35, were used. The lattice size was 
20^ X 64, or L « 2.6 fm. In addition, some prelim- 
inary results at a flner lattice spacing of a ~ 0.09 
fm are available. 

The big advantage of using lattices matched 
in lattice spacing and physical size for different 
numbers and masses of quarks is that the effects 
of dynamical quarks can be convincingly exposed. 
The simplest quantity that shows this is the static 
quark potential. Figure || shows the potential in 
quenched QCD and in three-flavor QCD with all 
the quarks at the strange quark mass. The slopes 
of these two potentials agree at the point chosen 
to define the length scale, but their overall shape 
is different. Notice that we do not see, and do 
not expect to see, string breaking in the distance 
range shown here. To quantify the change in 
shape we look at dimensionless quantities such as 
rQ^/a or ro/ri . Figure ^ shows ro for the three 
flavor runs. It also contains two flavor results with 
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Figure 6. Masses for 0++ propagators 

in quenched QCD (octagons) and full QCD 
(squares). The straight line is a linear fit to the 
large mass points, intended to represent the mass 
of a qq state, and the curved line is the mass of a 
two particle tt + rj state. 



Wilson quarks from CPPACS and SESAM, and 
a two-flavor improved Kogut-Susskind point from 
MILC. 

One aspect of lattice spectroscopy where sea 
quarks are expected to have important effects is 
in the chiral behavior of hadron masses. In the 
case of the Goldstone pion, at lowest order the 
squared pion mass is proportional to the quark 
mass, so corrections to this behavior can be dis- 
played by plotting the squared pion mass divided 
by the quark mass. Figure ^ is such a plot for 
the MILC runs in both full and quenched QCD. 
The increase at the right side is understood as 
the transition from chiral behavior to heavy quark 
behavior, while the much sharper upturn at low 
quark mass is the expected chiral logarithm. Al- 
though the quenched and three flavor curves look 
similar, the theoretical expectation is that the 
three flavor line has a flnite limit as ruq — *■ 0, while 
the quenched line diverges logarithmically. The 
quenched pion masses in Fig. ^ are similar to ear- 



lier quenched results of Kim and Ohta using the 
conventional Kogut-Susskind action at a lattice 
spacing of about 0.46 fm|lj]. A crude compari- 
son is in Fig. |^, where I have arbitrarily rescaled 
both the vertical and horizontal axes to match 
the results at the strange quark mass. 

While it is tempting to simply fit Fig. || to the 
chiral behavior predicted in Ref. [Q, it is prob- 
ably necessary to take account of the remaining 
flavor symmetry violations. Claude Bernard has 
been fltting this data both directly to the con- 
tinuum form, and to a form that takes into ac- 
count flavor symmetry violations using the La- 
grangian of Lee and Sharpe||T^ with empirically 
determined non-Goldstone pion masses in each 
loop, to calculate corrections to the pion mass[pO|. 
In the quenched case, where the coefficient of the 
chiral logarithm can be a parameter, reasonable 
fits can be obtained both with and without the 
fiavor symmetry breaking corrections. However, 
the coefficient of the logarithm comes out about 
0.06 with the simple fitting, and about 0.14 when 
the flavor symmetry breaking is included. When 
the three flavor data is flt, with the coefficients of 
the logarithms determined by the chiral theory, 
good fits can only be obtained when the flavor 
symmetry breaking is included. 

An application of these chiral corrections to 
the pion mass was presented at this conference 
by the Ohio State group pl|]. Here a particu- 
lar combination of chiral lagrangian parameters 
was computed and found to be inconsistent with 
niup = 0. The Boulder "HYP" action[|o) was 
used for the valence quarks in this project. In 
view of Bernard's observations on the effects of 
flavor symmetry breaking, it is appropriate that 
they used a valence quark action which makes 
these effects as small as possible. 

The MILC spectrum calculations on matched 
quenched and three flavor lattices found some 
interesting differences between the quenched and 
full QCD spectra (and at least one interesting 
lack of differences) . The most striking difference 
is in the coupling of hadrons to two particle in- 
termediate states, which is presumably present 
in full QCD calculations but is represented only 
by "hairpin" diagrams in quenched calculations. 
Figure || shows masses obtained for the isovector 
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Figure 7 
quarks. 

squares three flavor and the diamond a two flavor 
point. The burst is the real world value, and the 
cross the UKQCD quenched value [|2|. The hor- 
izontal scale parameterizes the quark mass, with 
lighter quarks to the right. 
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Figure 8. "Edinburgh plot" for MILC matched 
quenched and full spectrum calculations at a ~ 
0.13 fm. The squares are quenched, octagons 
three flavor, and the diamond two flavor results. 
The plusses are quenched points at a « 0.09 fm. 



0"*""*" propagator on matched quenched and full 
QCD lattices. We interpret the low quark mass 
behavior of the full QCD results as an avoided 
level crossing between a qq state and a two me- 
son Tr + rj state. In simple language, we are seeing 
a meson decay in the lattice simulations. 

The UKQCD collaboration pointed out that 
the ratio "J" , defined as 

dmv mK*im^-mp) 

J = a 2 ~ 2 2V' (1) 

drripg 2 (m|^ - m^j 

is a case where the quenched approximation 
produces an answer different from the real 
world [^. Figure shows this ratio on the 
matched quenched and full QCD MILC lattices, 
with a clear increase in the value when the sea 
quarks are included. This result can be com- 
pared with two flavor Wilson quark simulations 
by the CPPACS, UKQCD and JLQCD collabo- 
rations [^jl6[|2^, which also show a larger value 
in two flavor QCD. 



The nucleon to rho mass ratio is a much stud- 
ied quantity where lattice simulations typically 
disagree with the real world number. Many 
groups have looked at this, and it has become 
clear that this quantity is especially sensitive to 
effects of the lattice spacing. Thus it is inter- 
esting to use the matched lattices to look for 
effects of sea quarks. Figure || shows this ratio 
in the MILC calculations, with no discernable 
difference between the quenched and full QCD 
curves. This is in contrast to UKQCD calcu- 
lations on matched lattices with two flavors of 
Wilson quarks, where the two flavor numbers are 
larger than the quenched (See Fig. 9 of Ref. |p^ ). 
The two plusses in Fig. g are quenched points 
at a « 0.09 fm. Preliminary three-flavor re- 
sults at a w 0.20 and 0.09 fm show a similar 
trend. Although this lattice spacing dependence 
is much less than seen with the conventional 
Kogut-Susskind action, it is clear that a careful 
continuum extrapolation will still be required for 
this quantity. 



6 



Acknowledgements 

I am grateful to the organizers of Lattice-01 
for the opportunity to present this talk. I thank 
Claude Bernard, Joachim Hein, Takashi Kaneko, 
Francesco Knechtli, Kostas Orginos, and Howard 
Trottier for providing results used in this talk. 
This work was supported by the U.S. DOE. 



REFERENCES 

1. S. Naik, Nucl. Phys. B316, 238 (1989); C. 
Bernard et aZ.(MILC), Phys. Rev. D 58, 
(1998) 014503; C. Bernard et aZ.(MILC), 
Phys. Rev. D 55, 1133 (1997); G.P. Lep- 
age, Nucl. Phys. (Proc. Suppl.) 60A 267 

(1998) ; J.F. Lagae and D.K. Sinclair, Phys. 
Rev. D 59, (1999) 014511; Nucl. Phys. (Proc. 
Suppl.) 63, 892; K. Orginos and D. Tous- 
saint, Phys. Rev. D 59:014501 (1999); Nucl. 
Phys. B (Proc. Suppl.) 73, 909 (1999); K. 
Orginos, R.L. Sugar and D. Toussaint, Phys. 
Rev. D 60 (1999) 054503, Nucl. Phys. B(Proc. 
Suppl.) 83, 878 (2000).; Y.-B. Luo, Phys. 
Rev. D 57, 265 (1998); G.P. Lepage, Phys. 
Rev. D 59, (1999) 074502; C. Bernard et 
aL(MILC), Phys. Rev. D 61, 111502 (2000); 
U.M. Heller, F. Karsch and B. Sturm, Phys. 
Rev. D 60, (1999) 114502; F. Karsch, Nucl. 
Phys. B (Proc. Suppl.) 60A, 169 (1998); 
B. Beinlich, A. Bicker, F. Karsch, E. Laer- 
mann, A. Peikert, Nucl. Phys. (Proc. Suppl.) 
63 (1998) 895; O. Kaczmarek, F. Karsch, 
E. Laermann, Nucl. Phys. (Proc. Suppl.) 73 

(1999) 441; F. Karsch, E. Laermann, A. Peik- 
ert, Ch. Schmidt, S. Stickan, Nucl. Phys. 
(Proc. Suppl.) 94 (2001) 411. 

2. M. Golterman, Nucl. Phys. (Proc. Suppl.) 73 
(1999) 906. 

3. J. Hein, Q. Mason, G.P. Lepage and H. Trot- 
tier, these proceedings. 

4. P. Mackenzie, P. Lepage and H. Trottier, 
these proceedings; M. Nobes, these proceed- 
ings. 

5. C. Bernard et ai, Phys. Rev. D 61 (2000) 
111502; Phys. Rev. D 64 (2001) 054506. 

6. C. Bernard et al. Phys. Rev. Lett. 81, 3087 
(1998). 



Sara Collins, Robert G. Edwards, Urs M. 
Heller and John Sloan, Nucl. Phys. (Proc. 
Suppl.) 53, 877 (1997); private communica- 
tion; K.C. Bowler et al, Phys. Rev. D 62 
(2000) 054506. 



8. T. Blum et a/., |hcp-lat/0007038| ; Kostas 
Orginos, private communication. 

9. M. DiPierro and P. Mackenzie, these proceed- 
ings. 

10. A. Hasenfratz and F. Knechth, Phys. Rev. D 
63 (200 1) 114502; Phys Rev. D 64 (2001) 
034504; |hep-lat/0106014| ; F. Knechth, these 
proceedings. 

11. L. Levkova and T. Manke, these proceedings 

12. C. Bernard et al, Phys. Rev. D 64 (2001) 
054506. 

13. C. Bernard et al, Phys. Rev. D 62, 034503, 
2000. 

14. A. 



hep-lat/0105015; T. 



Ali Khan et al, 
Kaneko, private communication. 
15. G.S. Bah et al, Phys. Rev. D 62 (2000) 
054503. 



16. C.R. AUton et al, |hep-lat/0107021 . 

17. S. Kim and S. Ohta, Nucl. Phys. (Proc. 
Suppl.) 47 (1996) 350; Nucl. Phys. (Proc. 
Suppl.) 63 (1998) 185; Nucl. Phys. (Proc. 
Suppl.) 73 (1999) 195; Phys. Rev. D 61 
(2000) 074506. 

18. J. Gasser and H. Leutwyler, Nucl. Phys. 
B250 (1985) 465. 

19. W. Lee and S. Sharpe, Phys. Rev. D 60, 
114503 (1999) 

20. C. Bernard, in preparation 

21. G. Kilcup, G. Fleming and D. Nelson, these 
proceedings 

22. P. Lacock and C. Michael, Phys. Rev. D 52 
5213 (1995). 

23. JLQCD collaboration: S. Aoki et al, Nucl. 
Phys. (Proc. Suppl.) 94 (2001) 233; also these 
proceedings. 



